Flux-flow induced giant magnetoresistance in all-amorphous 
superconductor-ferromagnet hybrids 
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We present magnetoresistance measurements on all-amorphous ferromagnet (F) / superconductor 
(S) heterostructures. The F/S/F trilayers show large magnetoresistance peaks in a small field range 
around the coercive field of the F layers, at temperatures within and below the superconducting 
transition. This is interpreted as flux flow of weakly pinned vortices induced by the stray field of 
Bloch magnetic domains in the F layers. Bilayers show much smaller effects, implying that the 
Bloch walls of the F-layers in the trilayer line up and focus the stray fields. The data are used to 
discuss the expected minimum F-layer thickness needed to nucleate vortices. 



I. INTRODUCTION 



There are a number of phenomena currently under 
investigation which involve the combinations of super- 
conductors (S) and ferromagnets (F). In S/F/S config- 
urations, the superconductors can be coupled through 
the ferromagnetic layer, which may lead to so-called tt- 
j unctions i In F/S/F configurations, the superconducting 
transition temperature (T c ) depends on the relative ori- 
entation of the magnetization in two F layers^j^ which 
constitutes the so-called superconducting spin switch. In 
researching these phenomena, the question of the influ- 
ence of domain structures in the F layers is often ig- 
nored. Domain walls can have various effects. For in- 
stance, superconductivity can be enhanced by domain 
walls, through two different mechanisms. One is that 
Cooper pairs sample inhomogcncous exchange interac- 
tions in the wall, or the different directions of the mag- 
netization on the two sides of the wall, and experience less 
pair breaking. This was observed in bilayers of Nb and 
Permalloy (Py)£ and also in Nb/Coi The other mecha- 
nism comes about in ferromagnets with a preferred mag- 
netization direction perpendicular to the plane of a super- 
conducting film. Now, the presence of a domain wall can 
lead to a local reduction of the amount of flux through 
the superconductor, and therefore to less suppression (or 
relative enhancement) of superconductivity. This was 
demonstrated on S/F bilayers involving Nb and a fer- 
romagnetic garnet (BaFe^Oig)^, and on S/I/F bilayers 
(with I an insulating barrier) and F/S/F trilayers with 
Pb or Nb as the S layer combined with perpendicularly 
magnetized Co/Pd multilayers2iiS. 

A different situation occurs when the magnetization 
of the ferromagnet is in-plane and the magnetization in 
the domain wall rotates out of the plane (so-called Bloch 
walls). This can influence the superconductivity in the 
S layer if the flux from the wall creates vortices. Ob- 
servations on Nb/CuNi bilayers were interpreted in this 
waj*ii, but otherwise the problem has received little at- 
tention experimentally. Recently, the conditions for vor- 
tex formation were discussed theoretically (see ReiH and 



references therein) . In this work we present data from an 
experimental system well suited to observe the effects 
of vortices in transport measurements, consisting of a 
combination of an amorphous ferromagnet (a-GdigNisi, 
referred to as GdNi) and an amorphous superconductor 
(a-Mo2.7Ge, called MoGe). Because of the amorphous 
nature of the materials, the magnet has an extremely 
low switching or coercive field H c , corresponding to an 
applied flux density of less than 1 mT, while the super- 
conductor has very weak vortex pinning properties. Also, 
the magnetic material has a relatively high magnetiza- 
tion (due to the Gd atoms), which facilitates vortex for- 
mation. In F/S/F trilayers, we show the occurrence of 
extremely sharp resistance spikes when varying the mag- 
netic field around H c at temperatures near the base of the 
superconducting transition, which we interpret as due to 
flux flow. Furthermore, we find that the effect is much 
weaker in F/S bilayers, presumably because domain walls 
in both F-layers tend to line up, thereby focussing the 
flux coming out of the walls. Such a coupling of the do- 
main walls makes the F/S/F trilayer case different from 
the case of F/S bilayers. After the presentation of the 
data, we apply the model developed in Ref>iS to argue 
that vortices can be created in the MoGe layer by the flux 
from the domain walls, and we discuss the requirements 
of vortex formation in various other S/F systems. 



II. SINGLE FILM CHARACTERISTICS 

Our samples are grown on (100) oxidized Si by r.f. 
sputtering at room temperature, in a vacuum system 
with a base pressure < 2xl0 -6 mbar. Deposition rates 
were of the order of ~ 7.5 nm/min for the GdNi and 
~ 8.5 nm/min for the MoGe, as calibrated from low an- 
gle X-ray reflectivity. The compositions were found us- 
ing Rutherford Backscattering, and the amorphous na- 
ture of the films (i.e. the absence of crystallinity) was 
checked by X-ray diffraction. The bulk superconducting 
transition temperature T c of our MoGe films is about 
5.5 K, and such films show weak vortex pinning proper- 
ties as reported on similar material grown in the same 
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system 13 ' 14 . Another particular property of amorphous 
superconductors is that the very small mean free path 
(also reflected in a large specific resistance of typically 
about 200 x 10 _8 fhn) leads to a large zero-temperature 
London (magnetic) penetration depth Al(0), of order 
0.7 /mi. The zero-temperature coherence length £(0) of 
these films is small, around 5 nm. The numbers result in 
a quite small value for the zero-temperature lower criti- 
cal field H c i(0) of typically 1.3 x 10 3 A/m (correspond- 
ing to 1.6 mT), but in a very large value for the zero- 
temperature upper critical field H c2 (0) of ~ 13 T. 

Amorphous GdNi belongs to a general class of ferro- 
magnets combining a rare earth element and a transition 
metal element, which both carry a moment on their own 
subnetwork in the material. The amorphous state leads 
to a spatial distribution of the relative directions of the 
magnetic moments of both networks. If a net moment ex- 
ists, the state is called sperimagnetioi^. The moments of 
the two networks are coupled antiferromagnetically, and 
since the temperature dependence of the magnetization 
is different for both, there may exist a so-called compen- 
sation temperature T comp where the two magnetizations 
cancel. The case of Gdi-zNi^ is a special one. Accord- 
ing to the literatur e - 15 ; 16 , the Ni atom does not possess 
a magnetic moment below a critical concentration x c ~ 
0.8 while the Gd atoms have their full S-state moment of 
about 7 \iBi which order ferromagetically. With increas- 
ing x, for x < x c , the dilution effect of the non-magnetic 
Ni on the Gd matrix leads to a decrease of the magnetic 
ordering temperature T c down to 40 K at x — x c ~ 0.8. 
Above x c , a small moment appears on the Ni, antiparal- 
lel to the Gd moment, leading to a simple ferrimagnetic 
state. Obviously, a compensation temperature can only 
occur for x > x c . 
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FIG. 1: Magnetic moment M, normalised to its value at 4 K, 
vs. temperature T for a 20 nm GdNi film in an applied field 
of fioH — 3 mT. Inset: M vs. applied field ^lqH for the same 
film at T = 10 K. 

For our material with x = 0.81, we find a saturation 



magnetization, measured at 10 K on a film of a 75 nm, 
of 7.8 ± 0.2 x 10 5 A/m. Assuming the full Gd moment 
of 7 /is, this yields a small Ni moment of -0.02 [ib-, in 
good agreement with the earlier datait Fig shows the 
temperature dependence of the magnetization (normal- 
ized to the value at 4 K) of a 20 nm film, measured in an 
applied field of 3 mT. Clearly visible is that the magne- 
tization already develops around 220 K, but dips again 
to at T comp = 50 K. These value are reduced slightly 
for the thinner films used in the rest of this work, but the 
the alloy evidently orders at a much higher temperature 
than previously reported. It seems quite likely that T c 
was mistaken for T comp , possibly because of a too high 
applied field. The inset of Fig. ^ shows the field depen- 
dence of the magnetization, taken at 10 K. The coercive 
field H c is very small in this case, less than 0.1 mT, which 
is a consequence of the absence of anisotropy in the Gd 
S-state, and the lack of grain boundaries which hinder 
domain wall motion. 

All samples were patterned with e-beam lithography 
and broad beam Ar ion milling to 100 /im wide wires for 
a standard four point measurement with 1 mm between 
voltage contacts. The trilayer samples for which data 
is presented are of the form GdNi(x)/MoGe(y)/GdNi(z) 
with x, y, z all in nm, and the first of these grown on 
the SiO. All of the data presented are measured with the 
magnetic field applied in-plane and (anti-)parallel to the 
current to within a few degrees error (no precise align- 
ment procedure was undertaken), and with a constant 
current of ±100 fiA unless otherwise stated. 



III. RESULTS 

Fig. 13 shows the superconducting transition of a sam- 
ple GdNi(ll)/MoGe(21)/GdNi(ll), with a midpoint at 
2.16 K and a width (10 % - 90 %) of 60 mK. The re- 
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FIG. 2: Superconducting resistive transition for a patterned 
wire in a GdNi(ll)/MoGe(21)/GdNi(ll) trilayer. Inset: 
Magnetoresistance for H || J at 3.75 K. The arrows denote 
the direction of the field sweep. 
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duced T c indicates a significant proximity effect from the 
F layers. The inset shows the behavior of the resistance 
R versus applied field H for H | current / at 3.75 K 
(above the transition). Small dips are visible around the 
switching field of the F layers at 0.5 mT, which is the 
conventional anisotropic magnetoresistance with a mag- 
nitude AR/R max of the order of 9 x 10~ 5 . We observe 
a single peak in each quadrant, indicating that the two 
layers switch at the same field. Fig. |3| shows R(H) at 
2.15 K and at 2.095 K, at the base of the transition. 
Sweeping the field now leads to strong resistance peaks 
with AR several percent of the normal state resistance. 
For a related sample (with slightly lower H c ) we also plot 
the peak (dip) position H p d of the R{H) measurements 
through the superconducting transition (inset Fig. ■ It 
is clear that H c increases smoothly with T, and that the 
large magnetoresistance peaks in the transition are asso- 
ciated with the domain state of the F layers at H c . The 
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FIG. 3: Resistance vs. applied field of a 

GdNi(ll)/MoGe(21)/GdNi(ll) trilayer for two temper- 
atures within the superconducting transition as indicated. 
The arrows denote the direction of the field sweep. Inset: 
Values of the peak / dip field H p d in R(H) vs. temperature, 
for a related sample with slightly lower H c . 

maximum AR for this sample was — ' 26 r2, representing a 
change of 6 % of the total normal state resistance of the 
whole trilayer; however AR decreases with decreasing T 
and with this measurement current of 100 fxA the voltage 
in the peaks passes below a 1/zV criterion at 2.03 K. 

Since the resistance shows a peak in the the domain 
state of the F layers, we cannot interpret the data in 
terms of mechanisms which yield a decreased resistance 
(enhanced superconductivity) when the relative magne- 
tization directions in the two F layers deviate from par- 
allel. In principle, quasiparticle trapping could provide 
a mechanism for increased resistance^. It would then 
have to be argued that the domain state locally leads to 
antiparallel configurations between the two F layers, as 
was actually found in the case of F/S/F trilayers involv- 
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FIG. 4: Current-voltage characteristics at T = 1.8 K at an 
applied field —H c (drawn line) and slightly above +H C (dotted 
line) for the GdNi(ll)/MoGe(21)/GdNi(ll) sample. Inset: 
R(H) when biasing at 1 mA, above the depinning current. 
The arrows denote the direction of the field sweep. Also shown 
is a sketch of the sample configuration, with the directions of 
applied field, current, and flux inside the F-layer (including 
the domain wall) as indicated. 



ing (La,Ca)Mn03 and YP^CuaOy. 19 However, in view 
of the weak pinning properties of the superconductor, 
another possibility is flux flow resistivity associated with 
the motion of vortices formed spontaneously above Bloch 
domain walls. The sample is in the force-free configura- 
tion (H || /) for the applied field, but induced vortices 
pointing out of the plane of the sample will feel a Lorenz 
force due to the applied current, which can cause vortex 
motion across the width of the wire. To clarify this, we 
measured current (I) - voltage (V) characteristics on the 
same sample GdNi(ll)/MoGe(21)/GdNi(ll) at 1.8 K, 
well below the transition, and with the field either at - 
H c or sligtly above +H C . They are shown in FigQ] which 
also shows a sketch of the sample configuration, with the 
directions of applied field, current, and flux inside the 
F-layer (including the domain wall). At this tempera- 
ture we find a true supercurrent, and a gradual onset of 
voltage. Using a 1/iV criterion, the critical currents in 
the low and high resistive states taken from FigQ] are ~ 
640 and 340 respectively. Resistance peaks at this 
temperature therefore can still be seen, as long as the 
bias current is large enough to depin the vortices. This 
is shown in the inset of FigQJ where the voltage (resis- 
tance) was taken at a bias current of 1 mA. 

The effect is not strongly sensitive to variation of the S 
and/or F layer thicknesses. This is demonstrated by the 
data in Fig. 0in which MoGe layer is thicker, as well as 
the inset of Fig. [5] (thicker GdNi). These all show quali- 
tatively similar switching behavior to the original sample. 
In the case of the thicker MoGe layer some additional fea- 
tures are observed around zero field. These are the two F 
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layers switching independently, (also confirmed by mag- 
netoresistance measurements above T c - not shown here) . 
This is most likely a combination of a reduction of the 
direct coupling between the F layers for thicker MoGe, 
and an increase of the roughness and therefore H c of the 
top GdNi layer for a thicker spacer. With thicker GdNi 
the H c is reduced, the peaks can shifted to below 0.1 mT 
leading to a sensitivity in the switching at the steepest 
part of the curve (increasing field sweep) > 60il/mT at 
an applied field of ~ 90/iT. In Fig. [SJwe also show the ef- 
fect of patterning wires of different widths. The switching 
field changes from a 100//m to 2/Ltm width due to shape 
anisotropy, and with it therefore the field at which vor- 
tices are present in the MoGe. This allows us further 
control over H c in complement to varying the thickness 
of the GdNi. 




terface, and vortex pinning may actually be sensitive to 
(inhomogcncities at) the free interface. Also, given the 
observation of only one resistance peak in the trilayer 
case it appears that the domain walls couple across the 
MoGe layer, at least for relatively thin MoGe. We made 
similar observations of only one resistance peak for trilay- 
ers with two different thicknesses of the F-layer, such as 
GdNi(ll)/MoGe(21)/GdNi(22). We can speculate that 
this domain wall coupling both enhances the local flux 
density in the superconductor and sharpens the switch- 
ing behavior. 
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FIG. 6: R(H) for the bilayer MoGe(21)/GdNi(22) at 3.88 
K, and current directions as indicated; a positive current is 
defined as parallel to positive H. The arrows denote the di- 
rection of the field sweep. 



FIG. 5: R(H) curve for GdNi(ll)/MoGe(42)/GdNi(ll) 
trilayer at ~ 5K for two different wire widths. Inset: 
GdNi(38)/MoGe(21)/GdNi(38) switching below 0.1 mT at 
2.53 K. Arrows denote the direction of the field sweep and 
the use of the right-hand scale for the 2 /jm structure. 

A point of interest is that in bilayer samples the switch- 
ing behavior is less pronounced and also more compli- 
cated, since we now find asymmetry in the peak value of 
the voltage (resistance), but also asymmetry with respect 
to the current direction. Fields sweeps are shown in Fig. 
[S]for a bilayer sample MoGe(21)/GdNi(22), at a tempera- 
ture of 3.88 K, near the bottom of the transition. For pos- 
itive current (parallel to positive H) a jump rather than 
a peak is seen at +H C , and a rather broad peak occurs 
at —H c . We also observe an increasingly resistive back- 
ground (suppression of the superconductivity) at higher 
fields. For the other current direction the reverse is the 
case: a peak in R occurs at +H C when sweeping from 
negative H. Unraveling this behavior would need exten- 
sive study of the I-V characteristics, which will be for fu- 
ture work. We can, however, identify several differences 
between trilayers and bilayers. For instance, the bilayer 
is in fact asymmetric : the order parameter is strongly 
suppressed at the S/F interface, but not at the free in- 



IV. VORTEX FORMATION 

We have argued above that the magnetoresistive peaks 
we observe are due to vortices formed in the S layer above 
Bloch walls present in the F layer. To make this more 
compelling, we now estimate whether such vortices can 
actually be expected to form. For this we use the results 
from the model discussed recently by Burmistrov and 
Chtchelkatchev.— Basically, they find the amount of flux 
coming out of a Bloch wall of width S situated in a ferro- 
magnetic film of thickness dp characterized by a volume 
(saturation) magnetization M s , and from a free energy 
consideration calculate whether this can lead to the for- 
mation of a vortex in a superconducting layer of thick- 
ness ds, characterized by a London penetration depth 
Al, placed on top of the ferromagnet. For the case of 
a single domain wall they find, for given 5, dp, ds and 
Xl, the minimum or critical magnetization M c needed to 
overcome the lower critical field H c \ of the superconduc- 
tor. Since in our experiments M s is a materials constant, 
it is more convenient to write this condition in terms of 
a minimum thickness for the F-layer d™ m , which takes 
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the form 



XH. 



cl 



2X/S 

l-32G\/{ir 2 S) 



(n5) < 4A 
(irS) > 4A 



(1) 



Here, A = A 2 /ds if g?s < Al, otherwise A equals Al; H c \ 
is given by (<E>o/47T//oA 2 ) ln(A/£), with £ the Ginzburg- 
Landau coherence length and $0 the flux quantum; G « 
0.916 is the Catalan constant, and S 1 ./. units. For MoGe, 
the relevant values are Xl « 0.7 /<m, which for a 20 nm 
film yields A =24.5 /im; and £ « 5 nm, leading to i7 c i w 
1.8 A/m, an extremely low value which is due to the com- 
bination of a large bulk penetration depth and a small 
film thickness. For GdNi, the relevant characteristic are 
M s and S. As discussed above, M s = 7.8xl0 5 A/m (cor- 
responding to 0.98 T) is relatively large. Values for S 
are not exactly known, but the weak magnetocrystallinc 
anisotropy leads to large wall widths, which we take of 
order 1 /im. The relevant limit is then ttS <C 4A, and 
eq. n yields dPp m = 2.9 nm. For the thicknesses we use, 
and under the assumption of Bloch walls, the flux from 
a domain wall is therefore easily large enough to create 
vortices. 

This is the main conclusion from the calculation, but sev- 
eral more remarks are in order. First, it is interesting to 
note that, in this limit, dp m does not depend on A (apart 
from the logarithmic factor ln(«;)), since H c i oc 1/A 2 . 
Then, we have disregarded the effect of the in-plane ap- 
plied field. In terms of the model, this is allowed since 
ds << A, which means that the field fully penetrates 
without more than the vacuum contribution to the free 
energy. Experimentally, it can be noted that misalign- 
ment effects are apparently not relevant, since vortices 
are only created in significant amounts in the domain 
state of the ferromagnet. Making a rough estimate, an 
applied field of 2 mT (outside the flux flow peak) under 
a 1° misalignment yields an induction of 3 x 10 -2 A/m, 
again much smaller than the estimated H c \. Further- 
more, we note that, although the flux density from the 
domain wall is well above the lower critical field /xo-ffei, 
it is still much lower than the upper critical field [1qH c i- 
With a typical value of -fiodH C 2/ 'dT ss 2.5 T/K, and tak- 
ing T-T c w 100 mK, fi H c2 {T) is found to be w 0.25 T, 
very much larger than /j,qH c i. In summary therefore, 
the reason that vortices can be observed close to the 
resistive transition in our MoGe/GdNi system is that 
the amorphous ferromagnet combines a reasonably large 
magnetization with a large domain wall width, while the 
amorphous superconductor combines a large penetration 
depth with a high upper critical field. 

Given these different constraints, it is of interest to 
consider the possibility of vortex formation in various 
S/F systems which are currently under investigation as 
7r-junctions or spin switches, especially those based on 
Nb such as Nb/Py, Nb/Co or Nb/CuNi. The F layers 
in these combinations are qualitatively different, with 
Py having large M s and large (5, Co having large M s 
and small 5, and CuNi having small M s and larger 
S. Also considered can be YP^CuaGv (YBCO) and 



System 


MoM s [T] 


8 [/im] 


d^ m [nm] 


M)/(4A) 


a-MoGe/a-GdNi 


0.98 


1 


2.9 


0.03 


Nb/Py 


0.7 


1 


5.7 


16 


Nb/Co 


1.6 


0.05 


5.9 


0.8 


Nb/CuNi(50) 


0.1 


0.25 


19.1 


0.16 


Nb/CuNi(10) 


0.1 


0.5 


20 


1.6 


YBCO/LCMO 


0.75 


0.05 


51 


0.06 



TABLE I: Comparison of approximate critical thickness of F 
layer to achieve vortices above a Bloch domain wall for various 
S and F materials. The column entries are : the combination 
of materials, fioM 3 , the estimate for 8, the calculated d™ n 
and the relevant limit for using eq.Q Two cases are given for 
Nb/CuNi, one with the general S-layer thickness of 50 nm, 
and one with ds = 10 nm, as used in Ref^-. 



Lao.7Cao.3Mn03 (LCMO). For Nb, we use typical values 
of ds = 50 nm, Al = 50 nm, £ = 12 nm, for YBCO we use 
ds = 50 nm, A^ = 180 nm, £ = 2 nm. The different values 
for the ferromagnets are given in Table |TV| together with 
the computed value for d™ m . This is of course based on 
the assumption of Bloch and not Neel walls in such thin 
films, which may not be the case, but the numbers are 
instructive nonetheless. The table shows that the com- 
bination MoGe/GdNi actually yields the lowest value for 
d™ m due to the combination strong magnet / large do- 
main wall. Still, for the strong magnetsthe numbers do 
not vary overmuch, dp m is typically a few nm. For weak 
CuNi it is significantly larger, which is interesting in view 
of the observations of Ryazanov et alm^. They found flux 
flow behavior in the I — V curves of a Nb wire at in-plane 
applied fields around the coercive field of a block of 18 nm 
thick CuNi on top of a portion of the wire and ascribed 
the effect to vortices induced in the S layer due to Bloch 
domain walls in the CuNi. Taking into account that the 
Nb layer in their case was only 10 nm thick, the estimated 
value d r p m is 20 nm, which is roughly the thickness used 
in the experiment. The Bloch wall scenario for this ex- 
periment appears not unreasonable, since the prepared 
state is in-plane magnetized, while CuNi has a tendency 
to perpendicular magnetization in this thickness range, 
as found for CuNi/Cu multilayers 20 . The largest value 
for d™ m is found for YBCO/LCMO, which is due to the 
large value of ln(«) in this system. 



V. CONCLUSIONS 

We have demonstrated that amorphous F/S/F het- 
erostructures can show large magnetoresistance associ- 
ated with vortex motion in the S layer, induced by mag- 
netic domains in the F layers. This magnetoresistance 
can be several tens of Ohms change in a field step of a 
few tens of /xT due to the combination of weak domain 
wall and vortex pinning in these amorphous materials. 
We note that this effect can be a relatively simple test 
for the presence of Bloch domain walls in a ferromagnetic 
film. Also, the strong signals may provide a possibility to 
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combine magnetic domain and flux logic 21,22 in a flexi- 
ble way, since we have demonstrated that both the GdNi 
thickness and wire widths can be effective tools to tune 
the fields at which the peaks in flux flow resistivity are 
observed. That said, we also should mention some prob- 
lems open for further research. One point we have not 
touched is the obvious question whether the measured 
increase in resistance can be tied to flux-flow resistivity 
Pff in a quantitative way from the standard formula 
Pff = PnH/H C 2i with p n the normal-state resistance. 
At the moment we cannot answer that question since, 
apart from the fact that our measurements have not been 
performed in the linear regime of the I(V) characteris- 
tics where homogeneous flow can be assumed, we know 
neither the local field, nor the amount of vortices (de- 
termined by the domain wall width) between the voltage 
contacts. Also disregarded in the discussion are possi- 
ble geometrical effects which would lower the entry field 
for vortices due to an inhomogeneous current distribu- 
tion. This touches different questions such as whether 



this allows smaller values than d™ m , but also whether 
the nucleation of magnetic domains, through their cre- 
ation of vortices, actually is facilitated by edges or cor- 
ners. Fabrication of structures with artificial nucleation 
points would be an interesting extension of the present 
work. 
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